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Abstract 

Background: Monitoring joint angles through wearable systems enables human posture and gesture to be 
reconstructed as a support for physical rehabilitation both in clinics and at the patient's home. A new generation of 
wearable goniometers based on knitted piezoresistive fabric (KPF) technology is presented. 

Methods: KPF single-and double-layer devices were designed and characterized under stretching and bending to 
work as strain sensors and goniometers. The theoretical working principle and the derived electromechanical model, 
previously proved for carbon elastomer sensors, were generalized to KPF. The devices were used to correlate angles 
and piezoresistive fabric behaviour, to highlight the differences in terms of performance between the single layer and 
the double layer sensors. A fast calibration procedure is also proposed. 

Results: The proposed device was tested both in static and dynamic conditions in comparison with standard 
electrogoniometers and inertial measurement units respectively. KPF goniometer capabilities in angle detection were 
experimentally proved and a discussion of the device measurement errors of is provided. The paper concludes with 
an analysis of sensor accuracy and hysteresis reduction in particular configurations. 

Conclusions: Double layer KPF goniometers showed a promising performance in terms of angle measurements 
both in quasi-static and dynamic working mode for velocities typical of human movement. A further approach 
consisting of a combination of multiple sensors to increase accuracy via sensor fusion technique has been presented. 



Background 

Recently, a novel type of wearable sensor capable of 
detecting strain fields has been proposed [1-11]. Textile 
based deformation sensors can be produced by coating 
a thin layer of piezoresistive material on conventional 
fabrics [1,4,6] or by knitting conductive yarns with non- 
conductive yarns [2,3,7,9]. In other works, conductive 
threads are stitched [10,11] or attached [5,8] to the top 
of the fabric. The main features of textile deformation 
sensors are flexibility and the preservation of the mechan- 
ical properties of the garments on which they are applied. 
Textile deformation sensors have several advantages com- 
pared to solid-state sensors: negligible weight, thickness 
and possibility of spreading a high number of measuring 
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points over a flexible substrate. Sensing garments can be 
designed by applying sensor strips to specific locations 
on normal cloth. Changes in body shape and/or geometry 
due to human movements can in principle be estimated 
by reading variations on the measured strain. 

However, it may be difficult to recover the relation- 
ship between fabric strains and biomechanics parameters, 
such as tri-dimensional geometry or angles. Several solu- 
tions have been proposed to address this issue. In [12] a 
multivariate interpolation on a grid of sensors was used. 
Laviola in [13] reviewed the algorithms for hand pos- 
ture recognition; Gibbs and Asada [5] described a knee- 
sensing garment made with conductive fibers attached 
to flexible skin-tight fabrics. Mattmann et al. [14] com- 
bined a supervised learning algorithm with conductive 
thread sensors for the detection of torso movements. In 
[15] coated sensing fabrics are obtained by the integration 
of conductive elastomer (CE) materials on textile fab- 
rics by an ad hoc screen printing procedure with variable 
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topology. An application of CE sensors aimed at detect- 
ing the upper limb movement for neurological application 
is described in [16-18]. Generally, CE based sensing gar- 
ments perform well for slow and wide movements, while 
the accuracy, transient time and hysteresis has limited 
their use in reconstructing fast and small movements, 
such as anatomical torsion. In addition, retrieving sub- 
ject posture from textile strain measurements is highly 
affected by the relative position of sensors with respect to 
the joint being monitored. This issue has been addressed 
with the use of tight-fitting garments, which can reduce 
user acceptance especially in home rehabilitation con- 
texts. Even in the case of adherent garments, it is not 
easy to obtain reproducible results due to the inevitable 
sliding/bending of the sensors on the textile and to the 
difficulty of wearing the garment in the same way after 
donning and doffing. These latter aspects have limited 
joint angle tracking due to the necessity of using a com- 
plex and long lasting calibration procedure [6] or have 
restricted sensing garment usage in gesture classification 
applications [8]. 

In [19], CE sensors were configured in a double layer 
structure capable of direct angular measurements for 
application in rehabilitation and biomechanics. Since 
these double layer angular sensors are less sensitive to 
precise positioning and to their intermediate bending pro- 
file, they have the potential to solve some of the issues 
described for textile based strain transducers. 

In [7,9,20], Knitted Piezoresistive Fabric (KPF) sen- 
sors were demonstrated as a good tool for biomechanical 
and cardiopulmonary data acquisition and constitute an 
improvement on CEs. Compared to CE, KPF materials 
perform better in terms of response time, making them 
more suitable for wearable motion-capture applications. 
In fact, the transient time of CE sensors is very long and 
requires dedicated algorithms for predicting the final out- 
put after solicitation [6,21]. In addition, producing CE 
sensors entails using trichloroethylene with consequent 
rigid constraints on manufacturing sensors. In [22], a data 
glove based on KPF used as strain sensors was devel- 
oped and successfully tested in monitoring human hand 
gestures. 

In this work we exploited KPF technology to create 
a new generation of wearable goniometers inspired by 
the methodology introduced in [19]. The use of a new 
material led to a new theoretical approach compared to 
previous works aimed at improving the performance of 
CE sensors. 

Methods 

Single and double layer angular sensor working principle 

This section describes the basic theoretical aspects of 
textile-based angular sensors. Both single and double layer 
configurations were analyzed according to the guidelines 



described in [19] for conductive elastomer sensors and 
differences between them are highlighted. 

In [19] the working principle of a single layer CE 
angular sensor is described extensively. Single layer (SL) 
sensors are made of piezoresistive films (i.e. with a neg- 
ligible thickness compared to length and width) attached 
to textile or other flexible substrates. The following 
assumptions form the basis of the analysis: (f) the sens- 
ing layer is characterised by isovolumetric deformations, 
{if) the material resistivity p is constant, and (Hi) the fab- 
ric/flexible substrate is not extensible under bending (i.e. 
one side of the film has a constant length / during bend- 
ing). Figure 1 shows the CE film (in light gray) attached 
to the inextensible and flexible substrate represented by 
the thicker black line of length /. When the sample is 
flexed, the bottom side attached to the inextensible sub- 
strate maintains constant length /, while the sample height 
and length change according to its flexion. 

Under these assumptions, we demonstrated that the 
layer electrical resistance R$l of a specimen of length / 
parameterized in its arc length s is a function of the total 
curvature Aa, defined as the integral of the local curva- 
ture k(s) in the interval 5 e [0, /] [23] through the following 
relation: 

Rsl = 1-jt ~ ^ Aa + O ( sup k(sf ) 
ah d \se(o,l) J 

= l 2 ^-- P -Aa + 0[ sup k(sf) 

v 0 a ye (0,0 / 

where Vq is the specimen volume, d its width, Aa is the 
angle between the tangent planes to the sensor extrem- 
ities, ho is the initial thickness of the specimen, and 
O (sup (/c(s) 2 )) is a second-order infinitesimal function 
which tends to zero if the curvature k(s) -> 0. Note that 
k(s) assumes high values only when the sensor rapidly 
changes direction (as in cusps), which makes the function 
O (sup s€ ( 0 ,/) ^( 5 ) 2 ) negligible in detecting human body 
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Figure 1 Single layer CE sensor (light gray) attached to the 
inextensible substrate represented by the black line of length /. 
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shape. According to equation (1), the assumption of con- 
stant length guarantees the dependence of R$l only on its 
bending. As underlined in [24], the dependence of Rsl on 
Aa is also not affected by the particular bending profile. 
By neglecting the second order error in k(s), the angle Aa 
can be computed by: 



dl 2 d 

Aa = — R S l 

Vo P 



(2) 



Double layer (DL) angular sensors are based on the 
measurement of the resistance difference (ARdl = Rli — 
Rli) of two identical piezoresistive samples coupled via 
an insulating layer. This three-layer device is shown in 
Figure 2 where L\ and L<i are the piezoresistive layers 
with the same rest thicknesses ho. Rli and Rl2 are the 
electrical resistances of the two layers, respectively. Lq 
represents the insulating layer whose thickness is neg- 
ligible compared to ho. Lo can be both stretched and 
bent. 

In the case of device deformation without bending, as 
represented at the bottom of Figure 2, the electrical resis- 
tance difference between the two layers remains zero, 
since both the sensing layers undergo the same transfor- 
mation. When the device is flexed (as in Figure 3), [19] 
demonstrated that there is a linear dependence between 
ARdl and Aa: 



AR DL = Rli Rli 



2^Aa + 0( sup k(s) 3 
d \se(0,l) 



(3) 



This configuration proves that the error in angle esti- 
mation is reduced to a third order infinitesimal function 
with the maximum of the curvature k(s) [19]. While 
the hypothesis of iso-volumetric deformation and con- 
stant resistivity are still needed, the inextensibility can be 
removed since in (3) there is no dependence on the actual 



sample length. By neglecting the third-order error term, 
the angle Aa can be estimated by: 



Aa = ——ARdl 
2p 



(4) 



Sensor development 

Textile-based sensors were produced using knitted 
piezoresistive fabrics (KPF) which contain 75% electro- 
conductive yarn (Belltron®, produced by Kanebo Ltd) and 
25% Lycra®, manufactured as single jersey in a circular 
machine as described in previous works [7,25]. Circular 
electronic seamless-wear knitting machines by Santoni™ 
[26] were used to produce piezoresistive fabrics due to 
their capability to handle yarns with high elastic recov- 
ery. Flat knitting machines are more sophisticate in term 
of stitch selections, but less efficient in the handling of 
elastic yarns and in the production time. Moreover a 
Santoni machinery requires only the use of 4/8 spools of 
yarns compared with warp knitting machines. A conduc- 
tive bi-component fiber yarn based on polyamide loaded 
with carbon particles is used in combination with lycra 
to make this sensor. Piezoresistive fabric sensors change 
the electrical resistance according to the strain; the varia- 
tion in electrical properties is due to the modification of 
the interconnections geometry inside the fabric structure. 
Usually this property can be observed in stretchable fab- 
rics where the elongation of the fibers affects the flow of 
carrier inside the structure. When the conductivity of the 
yarn is due conductive particles as in bi-component fibers, 
the elongation of the yarn affects the charge transport 
mechanisms. The interconnections between the fibers 
and stitches are altered by the applied deformation. The 
elongation of the fabric modifies the distance between 
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Figure 2 Double layer angular sensor at rest (upper drawing) and stretched (lower drawing). When the sensor is stretched, the thickness and 
length change their values {ho -> h' Q and Iq -> Iq + A/). 
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Figure 3 Bent double layer sensor. L] and L 2 are the piezoresistive 
layers. Lq represents the insulating layer. 



stitches as well as the arrangement of the fibers in the yarn 
leading to a different geometry of interconnections. 

The samples were created according to the single and 
double layer design introduced in the previous section and 
hereinafter indicated as SL and DL sensors, respectively. 
The SL structure is made by coupling one rectangular KPF 
sample with an elastic fabric through a double-sided adhe- 
sive membrane. The DL structure is produced by adding 
another identical KPF layer to the back of the elastic fab- 
ric using the same adhesive membrane of the SL. Each 
piezoresistive layer has four semicircular pads specifically 
designed for sensor wiring. Figure 4 shows the electri- 
cal connections made of textile conductive yarns with 
PVC insulation (Bekinox®, produced by Bekaert). After 
removing the insulating material in the thread extremities, 
the connections were fixed using an ultrasonic welding 
device. 




Figure 4 Detail of one side of a double layer KPF sensor focusing 
on the connections between sensor and conductive wires. 



Figure 5 shows the structure of the DL sensor, where an 
insulating layer Lq is placed between the two piezoresis- 
tive layers L\ and L 2 . The topology of the connecting pads 
was specifically designed for a four-point measurement in 
order to minimize the effect of contact resistances. A con- 
stant current / is supplied through the external pads and 
the voltages V L \ = Vb x b 2 and Vn = Vb 4 b 3 (Figures 5 
and 6) between the internal pads are measured and can be 
related to resistances Rli and Rl2 (by dividing by/). As the 
curvature angle is related to the difference between resis- 
tances Rli and Rl2> a dedicated electronic was designed to 
continuously acquire ARdl- 

The custom-designed measurement setup is shown in 
Figure 6. A current generator supplies the two series of 
impedances with a constant current (/) and the acquisi- 
tion system is a high input impedance stage constituted 
by two instrumentation amplifiers (Al and A 2). The out- 
puts are amplified and the difference is measured in the 
following stage (A3). Both the resistance difference and 
the single layer resistance values are then acquired by a 
National Instruments acquisition board. The goniome- 
ters (Figure 4) that we tested have an overall rest length 
of 100 mm, which corresponds to the distance between 
the current carrying pads. The rest distances between the 
voltage sensing pads B\B 2 and B4B3 are 50 mm. The rest 
width and thickness are 10 mm and 0.5 mm, respectively. 

Double layer device calibration 

Equation (3) highlights the advantages of independence 
on the length variation of the device and the reduced 
errors in approximating the devices electrical resistance 
via a Taylor expansion. On the other hand, it only holds 
if the two piezoresistive layers have the same electrical 
properties. This property can be difficult to obtain from a 
practical point of view and the advantage of coupling two 
layers may be lost. To compensate for differences between 
the two piezoresistive layers, a calibration procedure can 
be applied. Considering the equation (1) for the layers 
L\ and truncated at the second order term, as in the 
following: 

R Ll =l 2 ^--^Aa (5) 
Vo a 

and 

Ri* = ( /r ) 2 ^7 + 37 A <* = tfo + l 2 ^qi ~ ~, Aaq 2 (6) 
Vq a Vo a 

where qo, q\ and q 2 are three corrective coefficients which 
take into account the geometrical and electrical diversity 
between the two layers due to the production process. 
qo, qi and q 2 are invariant with respect to the curvature 
and can be estimated by comparing the electromechani- 
cal characteristics of the two sensing specimens which will 
be introduced below. By multiplying equation (5) by q 2 , 
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Figure 5 Functional structure of a double layer KPF sensor. The central pads are for voltage acquisition, while the extremities are used for 
current supply and for the connection of the two single layer sensors. 



subtracting (6) and solving with respect to Aa and /, the 
following expressions are obtained: 



Aa 



qiRn - Rl2 + ffo 

§(tfl+*2) 



and 



Vb 



cfeRn + Rli - ffo 
§(<7i+<72) 



(7) 



(8) 



Experimental set-up 

KPF samples were characterised in quasi-static condi- 
tions for elongation and flexion both for SL and DL 



configurations through dedicated bench testing setups. 
In addition, DL dynamic performance in terms of flex- 
ion were preliminary evaluated on a human subject in 
comparison with inertial measurement units (IMUs). 

Quasi-static elongation test 

Elongation response was assessed using a custom- 
designed electro-dynamic testing system based on a lin- 
ear motor controlled by a PLC and able to apply strain 
cycles with controlled amplitude and velocity. Each layer 
resistance (Rli, Rli) and their difference (ARdl) were 
recorded as a function of the applied elongation. As shown 
in Figure 7(A), one extremity of the sample was attached 




Figure 6 KPF goniometer electrical equivalent and block diagram of the electronic acquisition system. The high input impedance stage is 
based on two instrumentation amplifiers (A1 and A2). A3 is a differential amplifier. The voltage output AV is proportional to the resistance 
difference and thus to the flexion angle. 
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Figure 7 Quasi-static test set-up. Figure 7 A shows the quasi-static elongation test set-up. The sample is attached to two clamps, one fixed (on the 
right) and the other one free to move (on the left). The slide has an 1 1 cm maximum range. The tested sample was subjected to a total deformation 
of 5 mm through 1 1 equivalent steps, both in elongation and shortening. Figure 7B shows the experimental setup for flexion characterization. The 
sensor is attached to a flexible substrate and is coupled to a commercial electrogoniometer. One extremity of the structure was clamped in a bench 
vice and the other one remained free to pivot around. 



to a fixed point while the other extremity was free to move. 
The sample was subjected to a total deformation of 5 mm 
through 10 equivalent steps, both in terms of elongation 
and shortening. Each step lasted one minute and the aver- 
age value of the recorded data was computed by the last 
30 seconds of the corresponding step. 

Quasi-static flexion test 

Quasi-static flexion characterization was carried out by 
relating Rn, Rl2 and ARdl with the output of a commer- 
cial electrogoniometer. Electrogoniometers are commonly 
used as a gold standard for angle measuring in biome- 
chanical applications [27]. The KPF sensor was attached 
to a flexible substrate composed of woven fiberglass 
cloth and epoxy resin (i.e. standard printed circuit board 
material). Then, a two-axis electro-goniometer SGI 10 by 
Biometrics (±2°C accuracy) was attached to the oppo- 
site side of the flexible substrate. One extremity of this 
structure was clamped in a bench vice and the other 
remained free to pivot around, Figure 7(B). Starting from 
0° the structure was bent to 90° through 13 steps. In each 
step, the sensor was held to rest for about 60 seconds 
and the average value of the recorded data was computed 
within the last 30 seconds. The test was performed both 
in flexion (from 0° to 90°) and extension (from 90° to 0°). 

Dynamic test 

For a preliminary evaluation of the DL sensor perfor- 
mance in dynamic conditions, a double layer KPF sensor 
was applied to a knee band and compared with the out- 
puts of two IMUs (MTw provided by XSENS [28]) placed 
on the thigh and on the calf, and used as a gold standard 
measurement instrumentation. In this case, a goniometer, 
longer than the one used in the quasi-static tests described 
above (40 cm, distance between the internal pads 30 cm), 
was used to entirely cover the knee joint. A representation 



of the set-up is shown in Figure 8. Several tests were 
performed by moving the knee in controlateral monopo- 
dalic standing at different velocities (slow, medium 
and fast). 

Data analysis 

Stretching and bending data, acquired on SL and DL sen- 
sors using the setups described in the previous section, 
were analysed in order to assess the performance of: 

• SL in stretching (i.e. Rsl variation with respect to the 
applied strain) 

• SL in bending and estimate the error by applying the 
relationships (1) truncate at the second order term in 
Aa 




Figure 8 Dynamic test set-up: double layer KPF goniometer 
applied to a band for the detection of knee flexion/extension; 
two IMUs are fixed to the thigh and on the calf in order to provide 
a reference measurement of the knee flexion/extension angle. 
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• DL in stretching and estimate the error by supposing 
that both layers have identical behaviour 

• DL in bending and estimate the error by applying the 
relationships (3) trunked at the third order term in 
Aa 

and in addition to: 

• Estimate the parameters qo, q\ and q 2 introduced in 
(6) to calibrate the double layer device 

• Preliminarily assess dynamic performance of DL in 
bending 

Hereinafter, to simplify the notation, R denotes the 
resistance of a single layer device or of one of the two 
layers (instead of Rsl)> while AR indicates the resistance 
difference between the two layers (instead of ARdl) 

Single layer resistance vs. stretching characteristics 

A single layer piezoresistive device was tested to estimate 
the linearity of the electromechanical characteristic and 
evaluate the errors introduced by hysteretic phenomena. 
Ten elongation/shortening cycles were performed with a 
total elongation of 5 mm by steps of 0.5 mm. The average 
values and the standard deviations of the electrical resis- 
tance were calculated (equations 9 and 10 respectively) for 
each imposed deformation d{ over the set of cycles. 

*ii = \H**i* (9) 

P 

In (9) and (10) p denotes the p th trial executed for a cer- 
tain deformation d{ and P the total number of trials. To 
roughly estimate the SL electromechanical properties, Rj. 
vs. di characteristic was approximated by a linear function 
and the deformation sensitivity (Sd) was computed as the 
slope of its linear approximation. 

Using the single layer device as a goniometer 

Relationship (1) truncated at the first order term relates 
the angle Aa to the electrical resistance R. If the length 
/ of one side of the device is maintained constant by an 
external constraint (e.g. the layer is glued on an inex- 
tensible film), the bending angle Aa can be estimated 
using relation (2). In order to validate (2), fifteen flex- 
ion (increasing angles) and extension (decreasing angles) 
cycles were applied. In each cycle the electrical resistance 
was acquired for 13 different angles Aat in the 0° - 90° 
range using the set-up described above. The average val- 
ues R/±ai and the standard deviations a^ a . of the electrical 
resistance were computed for each angular value: 




where k is k th trial executed for a certain angle Aai and K 
the total number of trials. 

R/xat trend was approximated by a linear regression in 
the least square sense and the linear regression slope, i.e. 
the angular sensitivity of the single layer sensor SAaSL> was 
computed. 

Using the double layer device as a goniometer 

The double layer device was tested in stretching to prove 
AR independence on elongation and in bending to ver- 
ify the relationship between angles and resistance differ- 
ences. Elongation/shortening and flexion/extension cycles 
were the same as described for the single layer sensor. 

Data collected in the elongation/shortening trials were 
used to calculate mean (13) and standard deviation (14) 
of AR as a function of the deformation applied, without 
using the calibration procedure described in the corre- 
sponding section. 

AR di = lj2 AR di,P ( 13 ) 

P 

OAR di = ^E^^-A^) 2 (14) 

Flexion/extension data were computed to determine the 
relation between the calibrated sensor output AR* and 
the angle Aa. The calibrated sensor output (AR*) can be 
retrieved from equation (7): 

AR* = qi R L1 - R L2 - q 0 . (15) 

Also in this case the mean (16) and the standard devia- 
tion (17) of AR* were calculated for each angular position, 
and the angular sensitivity of the double layer sensor 
SacxDL was computed as the slope of the linear approxima- 
tion of the AR*Aai trend. 

a^a„« = ^E a ^ ( 16 > 




Parameters identification and calibration 

Relation (7) can be used both for calibrating the dou- 
ble layer device and for identifying the system parame- 
ters which are still unknown. The parameters involved 
in the reconstruction of the angle starting from the two 
resistance values are the sensor width d, its resistivity p 
and the dissimilarities between the two layers, which are 
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expressed by qo, q\ and q 2 . These five variables can be 
summarized in a restricted set as: 



a\ = 
a 2 — 

&3 



p(qi+q2) 

d 

P(<?l+<?2) 

go d 



p{qi+q2) 

which transform the relationship (7) into: 

Aa = a\ Rli + a 2 Rl2 + ci 3 



(18) 



(19) 



In order to determine the vector (a\ a 2 a 3 ) J ', a set of 
experimental trials were performed by imposing different 
angles on the device in Figure 4. After executing m mea- 
surements (corresponding to angles Aa;, i = 0..m), the 
data collected were included in the linear system: 



Rli, Rl2, 1" 

Rn 2 Rl2 2 1 



Ri\ m RL2 m 1 
If matrix A is defined as 







Aa 2 




_Aa m _ 





a\ 

a 2 
a 3 



(20) 



A = 



Rlu Rl2, 1 

Rn 2 Rl2 2 1 



(21) 



Ri\ m Ri2 m 1. 
the parameter values can be computed by 



a\ 

Cl2 
d 3 



(AA T y 1 A T 



Aai 
Aa 2 

Aa m 



(22) 



The relationship (22) was computed by considering a 
complete increasing/decreasing cycle corresponding to 26 
angular measurements in the 0° - 90° range, performed 
with the set-up for quasi-static flexion characterization 
described in Figure 7b. The identification procedure we 
are describing holds for all the goniometers we tested. 
In the following, since numerical considerations are nec- 
essary to evaluate the adopted method, we consider the 
behaviour of a particular device. The values of the param- 
eters for the chosen goniometer are: 



ai = 13 10~ 3 a 2 = -3 1(T 4 a 3 



-58.7164 



(23) 



There is a notable difference between \a\\ and \a 2 \, 
despite the fact that they are coefficients of the resistance 
of the two layers L\ and L 2 whose geometrical proper- 
ties are very similar. This may be caused by the numerical 
instability of the method used. In order to verify this, the 
conditioning number of A T A was computed: 



= \\A T A\\ 2 \\(A T A)- l \\ 2 = 7.5169 10 13 (24) 

and indicates that the matrix A is very unstable. The 
poor conditioning of the matrix can be caused by small 
angles between different rows (considered as vectors) or 
large differences in the norm of different columns. In 
order to prevent the small angles problem, selecting those 
related to very different values can reduce the number 
of equations. The parameters were identified by recom- 
puting equation (22) on a set of three measurements 
corresponding to Aa = 0°, 45°, 90°: 

-3 1(T 4 a 3 



1.3 1(T 3 a 2 



-78.5450 



(25) 



fjL 3 (A) = 1.0792 10 14 . (26) 

The conditioning number gets worse, thus the second 
strategy, consisting in acting on the columns of A, is 
needed. 

Note that a 3 in A (a\ a 2 a 3 ) T has the same value in 
each row and is a function of qo through equation (18). 
As shown in the Results section, since qo is strictly related 
to the average value of AR for the unbent sensor, a 3 can 
be estimated separately. If the third coefficient is obtained 
independently by evaluating the average AR value in a 
short time from a static trial corresponding to Aa = 0° 
the system (20) is reduced to 



(27) 



Aai 




'Rli, 


R L 2, " 




a 3 


Aa 2 




Rli 2 


RL2 2 




a\ 




a 3 












_a 2 _ 


+ 




_ Aa m _ 




. R Ll m 


R L2 m _ 




_a 3 _ 



consequently: 



a 2 



= (B T B T y 1 B T 





Aai 




a 3 






Aa 2 




a 3 






_ Aa m _ 




_a 3 _ 





(28) 



The numerical stability of the method, computed in the 
previous three angles, improves notably, as proved by 

fi 3 (B) = \\B T B\\ 2 \\(B T By 1 \\ 2 = 67.5488 (29) 

and the computed values hold 

ai = 1.059 10- 3 a 2 = -1.067 HT 3 . (30) 

The same coefficients, computed by considering only 
two measurements (Aa = 0°, 90°), are 

ai = 1.065 10~ 3 a 2 = -1.063 HT 3 . (31) 

with a conditioning number 

fi 2 (B) = 45.9070 (32) 

Numerical output (30) and (31) are characterized by 
the same values (they differ in non-significant digits) with 
very similar conditioning numbers. 
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Dynamic test 

The device in Figure 8 was calibrated according to the pro- 
cedure described for the calculation of parameters a\ and 
^2- A zero- value was acquired when the knee was com- 
pletely extended. IMUs were calibrated in the standard 
way, to virtually obtain the same orientation of the related 
frames in the same joint position. The flexion-extension 
angle deriving from IMU outputs was computed by taking 
into account the components of the rotation matrix which 
describe the orientation of the IMU frame on the calf with 
respect to the frame of the IMU placed on the thigh. The 
knee flexion-extension measured with the KPF goniome- 
ter Og(t) was compared with the angle obtained by the 
two IMUs 6 iMU if)* I n addition, the two statistics of (33) 
and (34) were used to perform an inferential comparison 
between the measurement systems 

X = 0 IM u ~ 0~g = f 1 OjMuif) - O g (t)dt, (33) 

h — to J t 0 



cr = \\0imu - 0 g \\ 2 = —-—rj \ 1 (OiMuit) - 0 g (t)) 2 dt. 
t\ - to y j tQ 

(34) 

A set of ten trials on different subjects was performed 
in three different conditions: slow-speed, medium-speed 
and high-speed flexions. For each trial, the statistics 
(33) and (34) were computed. A Students t-test was 
then performed to prove that the measurement samples 



obtained by the two different systems belonged to the 
same population (Ho -hypothesis, i.e. X is a zero-mean 
random variable). 

Results 

Single layer resistance vs. stretching characteristics 

Figure 9 shows the relation between the SL electrical 
resistance and the deformation applied. The average R 
values and the corresponding standard deviations ctr 
are reported in Figure 10 as a function of the defor- 
mation applied d. The deformation sensitivity Sj is 
11950 £2 mm~ l and the maximum standard deviation 
OMax = niax^. GR d , holds 5603 Q for 2.8 mm, correspond- 
ing to an estimated standard deviation in terms of length 
of 0.4 mm which is comparable with the amplitude of the 
steps applied. 

Using the single layer device as a goniometer 

Figure 11 shows the angle/resistance characteristic of the 
SL device and the corresponding standard deviations. The 
angular sensitivity SacxSL holds 763 Q /° and the maximum 
standard deviation OMax — maxAa; holds 6405 Q 
for Aa = 37°, corresponding to an estimated angular 
standard deviation of 8.3°. 

Using the double layer device as a goniometer 

Figure 12 shows the relation between the resistance 
difference of the double layer device and the applied 
deformation. 
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Figure 9 SL sensor: R vs. applied deformation for both elongation (blue lines) and shortening (red lines). 
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Figure 1 0 SL sensor: average R vs. applied deformation for both increasing and decreasing deformation cycles. The vertical bars represent 
two standard deviation units in length. 



Note that the baseline of Figure 12 represents the dis- 
similarity between the behaviour of the two layers in 
terms of strain. As can be observed from the quasi- 
constant behaviour of the curve in Figure 12, AR is not 
dependent on / (i.e. q\ & 1). For a first approximation, 



the average value of the function plotted in Figure 12 
holds qo and can be analogically subtracted by the sig- 
nal deriving from one of the two layers or from their 
difference. Considering the acquisition diagram shown in 
Figure 6, qo can be analogically subtracted by measuring 




10 20 30 40 50 

[Degrees] 



70 80 90 



Figure 1 1 SL sensor: average R vs. angle for both flexion and extension cycles. The vertical bars represent two standard deviation units in 
length. 
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Figure 1 2 DL sensor: average AR vs. applied deformation for both increasing and decreasing deformation cycles. The vertical bars 
represent two standard deviation units in length. 



+ #2 Aa 



(35) 



the initial offset of the un-bended sensor, which is done by ^* _ ^ _ _ / 2 p 

regulating the offset of one of the three amplifiers (Al, A2 \ <i 
or A3). After offset compensation, the relation between 

the angle Aa and the resistance difference AR can be Relation (35) was experimentally verified to evaluate 

obtained starting from equation (7) and (15) by imposing the error in angular measurement. Figure 13 reports the 

qi = 1 data obtained in 15 trials for 13 different angles from 




40 50 

[Degrees] 

Figure 13 DL sensor: AR* vs. Aa for increasing (blue lines) and decreasing (red lines) angles. 
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Figure 14 DL sensor: average Aft* vs. angle for both flexion and extension cycles. The vertical bars represent two standard deviation units in 
length. 



0° to 90°, acquired both in increasing and in decreasing S^aDL is estimated as 955 Q /° and the maximum stan- 

measurements. dard deviation a Max was evaluated 5100 Q for Aa = 37°, 

Figure 14 shows average A/?* and the associated stan- corresponding to an angular error of 5.3°. 

dard deviation plotted as a function of the the angle A further analysis was performed starting from the 

Aa. The angular sensitivity of the double layer sensor observation of the different sensor behaviour between 
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Figure 1 5 DL sensor: average AR* vs. angle for flexion (red curve) and extension (blue curve) cycles. The vertical bars represent two 
standard deviation units in length. 
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flexion and extension (Figure 13). The mean and stan- 
dard deviations of AT?* were evaluated for increasing 
and decreasing angles separately (AR* Aa , o^r*. and 

A^*Aa/> a AR* Aa . respectively) and plotted in Figure 15. 



AR 



Acii 



-Yar" 



*I/D 

k 



I/D 
Cr AR* Aa . 



_ J/D 



AR 



;I/D 



'DV 



(36) 



(37) 



The maximum standard deviations of the two charac- 
teristics were computed: 



°Max = ™^ a AR* Aai = 2 ' 5KQ 



J Max 



= max(7 A m A = 1.4 KQ 



(38) 



(39) 



which corresponds to the estimated angular standard 
deviations of 2.4° and 1.48° respectively. 

Parameter identification and calibration 

Figure 16 shows the results of the described identifi- 
cation methodology for an increasing/decreasing cycle. 
The output of equation (19), evaluated with coeffi- 



cients (23) and (31), was compared with the Biometrics 
electrogoniometer. In terms of angular reconstruction, 
there was no significant difference between the two 
methods. 

Given this result, the parameters can be reduced fur- 
ther. The knowledge that a — a\ — —a2 and qo = as 
guarantees the correct functionality of the device. The two 
values required can be computed by calibrating the device 
in two points. In this way a rapid calibration procedure 
can be executed directly on the body after the goniometer 
has been integrated into the garment. After subtracting 
the initial value qo computed as the average value of 
the resistance difference for Aa = 0, it is sufficient to 
acquire the sensor output for a known angle (e.g. 90°) to 
compute a. 

Dynamic test 

Figures 17 and 18 show two example comparisons 
between the DL KPF goniometer and IMUs during con- 
tralateral monopodalic standing tasks. 

Both in the slow and fast knee flexion-extension, the 
double layer KPF goniometer performed well in angular 
measurements (maximum error of 5°) and managed to fol- 
low dynamic knee movements for compatible velocities 
with those of human movement. 

To reinforce this first analysis, the statistical approach 
described in the Methods section was followed. 
Fixing a significance level as rj = 0.05, all the trials led to 
the same results, i.e. that there are no biases between the 
two measurement systems as reported in Table 1. 
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Figure 17 Dynamic comparison between the calibrated KPF goniometer and the IMUs during a controlateral monopodalic standing task 
characterised by slow knee flexion-extension movements. The red line represents the KPF goniometer output. The blue line is the knee flexion 
extension angle obtained by the two IMUs. Maximum error is 4.2°. 



Discussion 

Double layer KPF goniometers performed better than 
single layer sensors in terms of quasi-static angle recon- 
struction (5.3° vs. 8.3°). 

Commercial electrogoniometers, such as those used in 
our quasi-static comparison with a declared accuracy of 



±2°, are widely used in ambulatory measurements of 
the joint range of motion (ROM) and movement fre- 
quency/velocity/acceleration for both clinical and occu- 
pational evaluations [29-31]. Studies on goniometer 
accuracy have shown errors of a few degrees, with great 
dependence on the sensor positioning and on the cross 
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Figure 18 Dynamic comparison between the calibrated KPF goniometer and the IMUs during a controlateral monopodalic standing task 
characterised by fast knee flexion-extension movements. The red line represents the KPF goniometer output. The blue line is the knee flexion 
extension angle obtained by the two IMUs. Maximum error is 5°. 
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Table 1 Dynamic comparison between the KPF 



goniometer and IMUs in knee flexion-extension tasks 



Activity 


X 


a 


t 


H 0 


Slow flexion 


0.06 


0.28 


-1.56 


VERIFIED 


Normal flexion 


-0.4 


0.33 


0.93 


VERIFIED 


Fast flexion 


0.08 


0.40 


0.20 


VERIFIED 



The KPF goniometer and IMUs angular measurements were compared in a set of 
different trials for different subjects. For each trial, statistics introduced in (33) 
and (34) were computed and a classical Student's t-test was performed 
(H 0 -ipotesis, i.e. X is a zero-mean random variable). 



talk between joints [32]. In [33] the evaluation of wrist 
ROM through electrogoniometers was performed with an 
error of between 2.2° and 6.2° over a flexion-extension 
range of 80°. Considering the widespread use of elec- 
trogoniometers, the reported errors are considered to 
be acceptable in clinical practice. The above reported 
KPF goniometer performances are slightly worse than 
those of commercial electrogoniometers, however they 
are still comparable with the ones accepted in goniometry 
applications. 

KPF goniometers have also been compared with IMUs 
within dynamic trials. IMU-based joint kinematic esti- 
mations, widely described in [34], have a reconstruction 
accuracy that is lower than 3° for flexion-extension joint 
movements [35,36], making the good agreement of our 
dynamic test very promising. 

The results reported in Figure 15 confirmed that the 
performance in extension (1.48°) is more accurate than 
the flexion performance (2.4°). It should also be noted 
that the angular error of the double layer sensor is 
mainly due to a hysteretic phenomenon pointed out in 
Figure 15 (i.e. maximum distance between the increas- 
ing and decreasing curve holds 10.4 KQ corresponding 
to 10.9°), since the standard deviations (38) and (39) are 
consistently smaller than those of the double layer sensor 

(&Max)> 

Given this finding, the following approach was con- 
ceived in order to improve the overall accuracy in terms 
of biomechanical reconstruction with wearable goniome- 
ters. As a practical example we will now consider the 
problem of hip flexion-extension detection, although the 
proposed method could be generalized to many joint posi- 
tion reconstruction problems. Two different goniometers 
can be placed on the anterior and posterior sides of the 
pelvis girdle (e.g. corresponding to the proximal parts 
of the rectus femoralis and the gluteus maximus mus- 
cles), as shown in Figure 19. The hip flexion/extension 
movement is acquired by measuring the mean value AR M 
between the anterior (AR A ) and posterior sensor {AR P ). 
This approach requires an on-body calibration, since the 
first angle between the abdominal region and the ante- 
rior part of the thigh and the second angle between 




Figure 19 Proposed configuration for hip flexion-extension 
detection. Two goniometers are placed on the anterior and posterior 
sides of the pelvis girdle to improve the accuracy in terms of 
biomechanical reconstruction. 



the lumbar region and the surface of the hamstrings are 
not equivalent due to body structure differences. In this 
case an on-body calibration is needed to ensure that 
the two goniometers perform equivalent measurements 
(e.g. set the 0 — value for each goniometer when the 
hip is completely extended). After this correction, the 
two sensor outputs are comparable and can be used to 
mutually integrate their information. In this approach, 
one of the two sensors will not return the effective angle 
between the two surfaces where it is applied, but the 
bias introduced will provide the re-calibrated sensor with 
the same behaviour as the other one in numerical terms. 
Thanks to this configuration, when one of the two sen- 
sors exhibits an increasing response, the other one will 
be decreasing. Our method reduces the error due to hys- 
teresis: the new AR M characteristic will be the mean of 
the increasing and decreasing responses and the over- 
all standard deviation will be reduced, in the worse case, 
to the sum of the increasing and decreasing standard 
deviations. Figure 20 reports the AR M characteristic, 
obtained by combining the increasing and decreasing 
responses. 

The angular sensitivity was estimated in 960 Q /° and is 
almost equal to that of a single goniometer. The maximum 
standard deviation is 3470 Q which corresponds to an 
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Figure 20 AR M vs. angle for both flexion and extension cycles in the advanced configuration based on two KPF goniometers for hip 
flexion-extension detection. The vertical bars represent two standard deviations in length. 



estimated angular standard deviation of 3.6° and repre- 
sents a consistent improvement. 

Conclusions 

We have developed a novel type of wearable goniome- 
ter based on KPF technology. In this paper, the working 
principle and the theoretical approach of single and dou- 
ble layer configurations have been described. On the basis 
of this theory, sensors were designed and produced on a 
fabric substrate. A calibration procedure that takes into 
account the dissimilarity between the two layers was pro- 
posed. Both single layer and double layer goniometers 
were tested in quasi-static conditions and compared with 
standard instrumentation. Double layer KPF goniometers 
performed better in terms of angle reconstruction com- 
pared to single layer ones (5.3° vs. 8.3° maximum error). 
In addition, double layer sensors are not sensitive to elon- 
gation and thus they are more suitable for applications 
in wearable motion detection. A preliminary dynamical 
evaluation showed how the quasi-static results could be 
extended in dynamic conditions. In addition, we demon- 
strated the overall sensor performance could be further 
improved through the fusion of two KPF goniometers 
per joint . Future work will focus on extending dynamic 
modeling and testing. 
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